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The ground state of a rotating Bose-Einstein condensate 
with attractive interaction in a quasi-one-dimensional torus is 
studied in terms of the ratio 7 of the mean-field interaction 
energy per particle to the single-particle energy-level spacing. 
The plateaus of quantized circulation are found to appear if 
and only if 7 < 1 with the lengths of the plateaus reduced 
due to hybridization of the condensate over different angular- 
momentum states. 
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The Hess-Fairbank effect [Q— d isappearance of the an- 
gular momentum (AM) of liquid helium 4 as it is cooled 
down to absolute zero with its container kept rotating 
slowly — is an analogue of the Meissner effect in super- 
conductivity, and it may therefore be regarded as a hall- 
mark of superfluidity. The essential requisites for the 
appearance of this effect are the single- valuedness of the 
wave function and the presence of a single Bose-Einstein 
condensate (BEC). Recent realization of BEC of lithium 
7 lH has opened up new possibilities associated with the 
attractive interaction between atoms; here the Pock ex- 
change interaction could energetically favor the forma- 
tion of hybrid BECs, which might modify the quantiza- 
tion of circulation and the Hess-Fairbank effect. In this 
Letter we investigate these possibilities in terms of the 
conceptually simple geometry of a quasi-one-dimensional 
torus. 

We consider a system of N weakly interacting bosons 
with mass M, confined in an torus of radius R and cross- 
sectional area S — 7rr^, where for simplicity we assume 
r <^ R. This condition justifies our assumption that 
the radial wave function is fixed and independent of uj 
— the angular frequency of rotation of the torus. At suf- 
ficiently low temperature, the interaction between dilute 
hard-core bosons is well approximated by Fermi's contact 
interaction, which is characterized by the s-wave scatter- 
ing length a. The associated mean-field interaction en- 
ergy per particle is given by gN, where g = 2ah'^/MRS. 
The positive (negative) sign of g implies that the effective 
interaction between bosons is repulsive (attractive). The 
Hamiltonian of our system in the rotating frame is given, 
up to terms which are constant in our approximation, by 
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where Uc = h,/2MR^ is the critical angular frequency, 
Z,m, and n denote the projected angular momenta in 
units of h, and c| and c; are the creation and annihi- 
lation operators of bosons with AM /. In Eq. (|l|), we 
have added the term ?iWc(<x'/2wc)^c/c| = Nhuj"^ / [Aojc) 
which is compensated for by the Lagrange multiplier a in 
Eq. (^) and therefore does not modify any result below. 

We determine the minimum-energy state of the Hamil- 
tonian (|]) within a Hilbert subspace given by |^')hf — 
I • • • , n_;, • • • , n_i, no, ni, • • • , n/, • • •), where n; denotes 
the number of bosons that occupy the state with AM /. 
This is nothing but the Hartree-Fock (HF) approxima- 
tion; other possibilities will be discussed later. Because 
the total number of bosons is A'', nj's should satisfy 
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The expectation value of the Hamiltonian with respect 
to the state |4')hf is given by 
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where Ki{uj) = hujc{l ~ w/2wc)^- The distribution of 
{ni} is determined so as to minimize E{{ni}) subject to 
condition (||). 

Case of repulsive interaction. — We first show that our 
ansatz wave function |5')hf reproduces some well-known 
results. When 5 > 0, it is possible to simultaneously 
minimize the kinetic energy and the interaction energy 
in Eq. (^) with ni — N ii I = [{uj + u!c)/2lUc] and ui — 
otherwise, where the symbol [x] denotes the maximum 
integer that does not exceed x. This result implies that 
a single BEC is energetically favorable whether or not 
it is rotated; Bogoliubov's virtual-pair excitations only 
cause a depletion of the condensate and do not alter this 
conclusion. The single valuedness of the wave function 
dictates that the projected AM be quantized in units of h, 
but one needs something more to show that it is quan- 
tized in units of Nh. The Onsager-Feynman condition 
for the quantization of circulation, in fact, requires the 
more stringent latter condition. For the case of repulsive 
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interaction, the Fock exchange interaction favors a sin- 
gle BEC Q , thereby enforcing sharp transitions between 
different AM states and requiring that the circulation be 
quantized in a uniform system as considered in this Let- 
ter. (In a related context, Castin and Dum have recently 
considered the stability of vortices for the nonuniform 
case of parabolic potentials See also Refs. 

Case of attractive interaction. — When g < 0, it is 
impossible to simultaneously minimize the kinetic en- 
ergy and the interaction energy; the kinetic energy be- 
comes minimal when the distribution {ni\ peaks sharply 
at I = [{ijj + ti)c)/2ti'c], whereas the interaction energy 
becomes maximal for this case. Were it not for the ki- 
netic term, the lowest-energy state would be the one in 
which the distribution of n; is maximally spread; no sin- 
gle state I could then be macroscopically occupied, and 
there would be no BEC. When the system is spatially 
confined, however, the kinetic term competes with the 
attractive interaction, allowing a metastable condensate 
to be formed. 

The minimal-energy distribution {ni\ is determined so 
as to minimize E({ni\) in Eq. (|^), subject to condition 
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where a is a Lagrange multiplier, and 7 = \g\N / {Tiojc) — 
4:N\a\R/S is the ratio of the mean-field interaction en- 
ergy per particle to the single-particle energy-level spac- 
ing. To find an estimate of 7, we rewrite it as 7 ~ 
4 X l{)^'^N\a\[A]R[^m]/ S[^lm'^]. For the case of lithium 
7 with \a\ ~ 14. 6A, R = Ifim, and r = 0.2/iTO, we have 
7 ^ 0.046A^. With suitable choice of these parameters, 
it is possible to prepare the system both with 7 < 1 and 
with 7 > 1. 

For Hi to be positive, there must be minimum and max- 
imum values of I, i.e., —li and 12- Equation then be- 
comes = which upon substitution of Eq. (^) 
for ni gives 
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where uj = Lojlio^. With the definitions of l\ and ^2, we 
have (/i + < a < (/i + 1 + o))^ and (^2 - tD)^ < a < 
(^2 + 1 — ci))^, which lead to 

[ 4/1-2/2-1 , , 4^2-2/1-1 \ 7 
(/i -1-/2) max< V Lo^ LO} < 



/1 + /2 + I 



^n^i^o^ ■ [ 4^1-2/2 + 3 , , 4/2-2/1 + 3 
< (/1+/2 + 2) mm< h cj, uj> (6) 
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These inequalities uniquely determine the pair of integers 
(/i, I2) for a given set of 7 and uj. 

When the torus is at rest (i.e., w = 0), Eq. (^) becomes 
ni — {N/j){a — P), where a is given from Eq. (H) with 



/i = /2 as a = 7/(2/1 -f 1) + /i(/i + l)/3, and Eq. (|) 
reduces to /i(4/? - l)/3 < 7 < (/i l)[4(/i + if - l]/3. 
These inequalities uniquely determine the number 2/i + 1 
of macroscopically occupied AM states for a given 7. For 
example, 7<1,1<7<10 and 10 < 7 < 35 correspond 
to 2/1 + 1 = 1, 3, and 5, respectively. Thus, there is 
a single BEC when 7 < 1. This condition agrees with 
the usual criterion for a metastable BEC to exist that is 
obtained for a parabolically confining potential using the 
Gross-Pitaevskii (CP) equation [Q. A new finding in our 
analysis is that for 7 > 1 BEC becomes hybridized over 
different AM states. 

At the continuum limit Wc (i.e., R — > 00) with 
\g\N held constant, 7 and /i become infinite with a/7 ^ 
0(l//i). It follows from the relation rt; = {N/^){a — 
that all njs becomes vanishingly small, of the order 
of N/li. Thus, no BEC exists for an infinite system in 
accordance with the standard widsom ||^ . 

The analysis for the case of w ^ is straightforward, 
and we describe here only the results that are relevant to 
later discussions. 

1. The region in which a single BEC with ui — N exists 
is given from Eq. (^ with /i = — /, /2 = / by 
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When uj/lUc is an odd integer, condition (^ can never be 
met, so no unique BEC can exist no matter how weak 
the attractive interaction. 

2. The region in which two states with AM / and / + 1 
are macroscopically occupied is given from Eq. (|^) with 
/i = -/,/2 = / -t- 1 by 

I — -2/-l|<7<min|3 — -6/ + l,-3— + 6/ + 71 , (8) 
and the corresponding distribution of bosons is given by 
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3. In general, the region in which k states with AM 
/,/ + l,---,/ + fc— 1 are macroscopically occupied is given 
from Eq. (||) with /i = — /, /2 = ^ + fc — 1 by 

(-6l-2k+l UJ 6l + 4:k-5 UJ ] 
. r-6/-2fc+5 UJ Ql + Ak-l UJ ^ , , 

The phase diagram is shown in Fig. |l|. We have thus 
shown that there are regions of 7 and uj in which more 
than one AM state is macroscopically occupied. This pre- 
diction can be tested most directly by switching off the 
trap potential and let the system expand. Due to Heisen- 
berg's uncertainty relation, the tight radial confinement 
of the trap causes the gas to expand more rapidly in that 
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direction than in other ones, and the superposition of 
BECs having different AM should result in an interfer- 
ence pattern with broken axisymmetry ||] . 

Partial quantization of circulation. — When we fix 
7 = \g\N/?iu;c < 1 and increase uj from 0, we alterna- 
tively pass regions in which one or two AM states are 
macroscopically occupied (see Fig |^), and in the latter 
regions the distribution of bosons between the two AM 
states changes continuously with w, as can be seen from 
Eq. (|^). What happens then to the circulation k of the 
system? When a single BEC with AM I exists, k is given 
by hl/M. When two AM states are macroscopically oc- 
cupied, K should be given by h{l)/M, where (1) is the 
ensemble-averaged value of the AM. To find this value, 
let us restrict ourselves to the region 



{u!-UJc)/uJc\ < 7 < -3|(CJ - LUc)/llIc\+4:, 
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where two states with AM I — and I — 1 are macro- 
scopically occupied, and the number of bosons in each 
condensate is given from Eq. (^ by no = N/2 — N{uj — 
ujc)/{2ujcl) and ni = N — uq. Hence the ensemble- 
averaged AM (l) is given by 
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which does not show any sharp transition (see Fig. |^), 
in sharp contrast with the case of repulsive interaction. 
Suppose now that we perform the Hess-Fairbank exper- 
iment for the frequency of rotation and for the strength 
of interaction that satisfy the condition (|ll|). Then the 
AM will not completely be "expelled" even at absolute 
zero and have a nonzero value given by Eq. (p^). Only 
when those parameters are in the region (|) with / = 0, 
the AM should vanish at absolute zero. 

Hybrid BECs vs. a phase- coherent single BEC. — 
We have shown within the HF approximation that hy- 
brid BECs exist for some ranges of parameters 7 and 
cj/wc- Recently, Rokhsar has argued that hybrid or "frag- 
mented" BECs are inherently unstable against the for- 
mation of a single BEC whose macroscopically occupied 
state is a linear combination of the "fragments" with defi- 
nite relative phases |^ . In our situation, the "fragments" 
refer to macroscopically occupied AM states. We first 
discuss the stability of binary BECs against forming such 
a phase-coherent single BEC. Because the properties of 
the system are periodic functions of to with periodicity 
2ujc, we may consider, without loss of generality, the re- 
gion ( 11 ) in which two BECs with I = and I = 1 coexist. 
The state vector of this binary BECs is given by 
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where uq and ui are the numbers of bosons in the I ~ 
and I ~ 1 states, which are given by Eq. d). To 



be compared with this state is a single macroscopically 
occupied state whose creation operator 6^ is given by 
b^ = acjj + (3c\, where a and (3 are determined so as 
to minimize the total energy, subject to jap -I- = 1. 
The corersponding single BEC is given by 
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The crucial observation here is that when only two 
states are macroscopically occupied, the expectation 
value (J?) single of the Hamiltonian (Q) over the state ( p^ ) 
does not contain any non-HF terms that are of the same 
order of magnitude as the HF terms because of the con- 
servation of AM. Therefore, the system cannot lower its 
energy by establishing a relative phase coherence between 
the different AM states. The minimum value of (if)singie 
is reached when 
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and by a straightforward calculation, we find that 



H)iiF - (H) 
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However, because the energy difference is of the order of 
1/A^, the two states ( p^ ) and ( p4[ ) are virtually degen- 
erate. In real life, however, there are inhomogeneities in 
the container "walls" etc., which break the exact axisym- 
metry. Such a perturbation, however weak, could stabi- 
lize the single coherent BEC relative to a "fragmented" 
one. To show this, consider a symmetry-breaking per- 
turbation that mixes the I = state and the I — 1 
state: V = tc\ci + t*c\cQ. It is easy to see that while 
V does not lower the energy of the system for the HF 
state ((F)hf = 0), it does for the single coherent BEC; 
{V)^;ng\o = 2NRe{a*pt) = -2N\a*f3t\, provided that 
arga-arg/3-argt — ±7r. Because both I — and I = 1 
states are macroscopically occupied, i.e. a ~ 0(1) and 
P ~ 0(1), (V^)singiG is extensive. The single coherent 
BEC can therefore become energetically favorable due to 
a (possibly infinitesimal) symmetry-breaking perturba- 
tion. It should be noted, however, that the plot of (/) 
versus cu in Fig. ^ remains basically unaltered because it 
does not depend on whether or not a phase coherence 
is established between two macroscopically occupied AM 
state. 

The situation is different when more than two AM 
states are macroscopically occupied. Now the expecta- 
tion value of the Hamiltonian contains non-HF terms that 
arc of the same order of magnitude as the HF terms, so 
that without the need of the symmetry-breaking pertur- 
bation the system can lower its energy by establishing 
a relative phase coherence. To show this, let us con- 
sider the case of 1 -I- 3|w/wc| < 7 < 10 — 6|w/wc|, where 
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three AM states I = —1,0,1 are macroscopically occu- 
pied. The hybrid BEC state is described by |4')hf = 
|n_i, no, rii) and the corresponding single coherent BEC 
is described by |4')singie = (acLi +/3cJ + 7c|)^/ViVT|vac) 
with |ap + + I7P = 1. The expectation value of the 
Hamiltonian with respect to |4')hf is given by 

{H)uF = K-in-i + Kom + Kifii 

-|.g|(n_ino + n^ni + ni7i_i)- \g\N{N -l)/2, (17) 

which is minimized when nzpi — Af[l=p(3tiJ±Ci'c)/(7^c)]/3 
and no = ^^[1 + 2/7] /3. The expectation value of the 
Hamihonian with respect to |^')singic is given by 

(i?>si„gic - iV(i^_i|a|2 + + X1I7I') 

-|5|7V(7V-l)(|an/3|2 + |/3|2|^|2 + |^|2|^|2) 

-\g\N{N ~ l)/2 - \g\NiN - l){aP*'-f + a*(3^Y)- (18) 

Because the last two terms are phase-dependent and of 
the same order of magnitude as the remaining terms, 
it is clear that the single coherent BEC can have a 
lower energy than the fragmented BEC state by, e.g., 
the following choice of amplitudes, a — ^n_ie'^° , P = 
y^noe'^'', 7 = Y^nTe*^^ with the relative phase relation 
9a — '29f3 + 0^ = 0. Thus the ternary BEC state is un- 
stable against the formation of a single coherent BEC. 
Similar mechanisms should work when more than three 
AM states are macroscopically occupied. 

To summarize, we have studied the ground-state prop- 
erties of a rotating BEC with attractive interaction con- 
fined in a quasi-one-dimensional torus. When the con- 
dition (0) is met, only one AM state is macroscopically 
occupied. When the condition (||) is met, two BECs with 
different AM can, in principle, coexist. However, any de- 
viation from the exact axisymmetry is shown to stabilize 
a single coherent BEC relative to a "fragmented" one. 
The plateaus of quantized circulation appear if 7 < 1, 
but the lengths of the plateaus are reduced. In other re- 
gions of parameters 7 and w, more than two AM states 
are macroscopically occupied, where non-HF terms sta- 
bilize a single coherent BEC even in the presence of the 
exact axisymmetry. 
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